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Abstract. Quantum systems often contain negative energy densities. In general 
relativity, negative energies lead to time advancement, rather than the usual time delay. 
As a result, some Casimir systems appear to violate energy conditions that would 
protect against exotic phenomena such as closed timelike curves and superluminal 
travel. However, when one examines a variety of Casimir systems using self-consistent 
approximations in quantum field theory, one finds that a particular energy condition is 
still obeyed, which rules out exotic phenomena. I will discuss the methods and results 
of these calculations in detail and speculate on their potential implications in general 
relativity. 
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1. Introduction 

General relativity allows spacetimes of any geometry. Given any g^ u , we can compute 
and R, and then set up a matter configuration whose stress-energy tensor is 



to obtain a solution to Einstein's equations with the desired geometry. As a result, 
nothing seems to prohibit the existence of exotic phenomena such as closed timelike 
curves [1], traversable wormholes [2], or superluminal travel [3]. Therefore we expect 
there to exist some restrictions the possible tensors T^ u , usually called energy conditions, 
which would lead to restrictions on the possible spacetime geometries. While it is 
straightforward to show that classical field theories obey energy conditions that are 
strong enough to forbid exotic phenomena, quantum field theories appear to violate 
these conditions. 

We will focus on the following energy conditions: 

• The Weak Energy Condition (WEC) requires that for all timelike vectors V^ 1 , 



That is, all observers see positive energy density. 

• The Null Energy Condition (NEC) is weaker than the WEC, and requires that 
Eq. (2) hold only for null vectors V M '. 

• The Averaged Null Energy Condition (ANEC) is weaker than the NEC, and only 
requires that the NEC hold when integrated over a complete null geodesic. 

In a classical background, all of these conditions can be imposed consistently, and 
any one of them would be sufficient to rule out exotic phenomena [1] . On the other hand, 
the standard Casimir system of parallel conducting plates has a static negative energy 
between the plates, and thus it appears to violate all of these conditions. However, there 
are some obvious caveats to this result. An external agent has to hold the plates apart, 
and impose the boundary condition. Furthermore, in the case of ANEC, the geodesic 
has to go through the boundary, where we are ignoring the effects of the material. Since 
gravity couples to all sources of stress-energy, we cannot consider the contribution of 
the Casimir energy without also including the contributions due to the materials and 
their internal interactions. 

Previous work [4, 5] avoided these problems by considering a domain wall 
background that is renormalized with standard counterterms. In this way, all the 
contributions to the stress-energy tensor can be included. The result, shown in Figure 1, 
is that for any coupling, there is always a small region at large enough x where the 
energy density is negative. So we can always choose the coupling small enough that our 
semiclassical approximation is reliable and still find negative energies. This situation 
also violates quantum inequalities [6, 7, 8]. Although this system violates WEC, since 
we have a region of negative energy density, and NEC, if we take a null vector in the 
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Figure 1. Classical energy density due to the wall (dotted), the quantum correction 
(dashed), and the total (solid), for typical values of the coupling constants.. 



region of negative energy pointing perpendicular to the wall, ANEC is still obeyed, since 
the contribution from the complete geodesic perpendicular to the wall is only negative 
if the total energy of the domain wall is negative (indicating an unstable vacuum and a 
breakdown of our approximation). 

In a number of other examples, explicit calculation shows that ANEC is obeyed. 
These include a geodesic outside a spherically symmetric background potential [9] or 
a dielectric sphere [10]. Other calculations also show that energy condition violation 
is more difficult to achieve in realistic situations than idealized models would suggest 
[11, 12]. ANEC is also known to be obeyed by free scalar [13] and electromagnetic 
[14] fields in flat spacetime. Other works have found restrictions on energy condition 
violation in flat space [15, 6, 16]. 

2. Plate with a hole calculation 

To further investigate ANEC, work done in collaboration with Ken Olum [17] considers 
another alternative. We avoid the plate by drilling a hole in it for the geodesic to travel 
through. Then we might expect that the region around the hole would provide only 
a small correction to the negative contribution from the rest of the geodesic, yielding 
a violation of ANEC. Since the geodesic never encounters the material, the result will 
be finite with no contributions from the counterterms (which have support only where 
there is a potential). Our approach will be to use a Babinet's principle argument to 
show that the Casimir energy of a Dirichlet plate with a hole is the sum of the Casimir 
energy of a full Dirichlet plate and the Casimir energy of a Neumann disk. Then we 
will use scattering theory in ellipsoidal coordinates to solve the disk problem, and add 
in the standard result for the full plate. 

In free space, we can decompose the spectrum into modes that are odd or even 
under reflection of the z-axis. Imposing a Dirichlet boundary condition at z = has 
no effect on the odd modes, since they already obey the boundary condition. However, 
the even modes are modified: They turn into the odd modes multiplied by the sign of 
z, with a cusp at z — 0, as shown in Figure 2. 

If there is a hole in the plate, we find even functions that obey Neumann conditions 
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Figure 2. In free space (left) we have even and odd normal mode wave functions. 
In the presence of a Dirichlet plate (right), the even functions are replaced by the odd 
functions with a change of sign crossing the plate. 



in the hole (since they are even) and Dirichlet conditions elsewhere. Up to a similar 
sign flip between positive and negative z, these are the same functions as one obtains 
for the odd modes of a Neumann patch with the same shape as the hole, as shown in 
Figure 3. 




Figure 3. If there is a hole in the Dirichlet plate (left), the new even functions satisfy 
Neumann conditions in the hole and Dirichlet conditions elsewhere. The odd functions 
for a Neumann patch with the same shape as the hole (right) are the same except for 
a change of sign between sides. 



Let E be the contribution to the Casimir energy from the free even modes, O be the 
contribution to the Casimir energy from the free odd modes, and A be the contribution 
to the Casimir energy from the modes for the plate with the hole. Then we have the 
following total contributions to the energy, where we have taken the difference with the 
free space result: 

Free space: E + O subtract jree space miQ 

Dirichlet plate: O + O subtract jree space Q _ E 

Dirichlet plate with hole: A + O subtract jiree space A _ E 

Neumann disk: E + A subtract free space A _ Q 
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Thus we have found 

[Dirichlet plate with a hole] = [Neumann disk] + [Dirichlet plate] . (3) 

Next, we set up the disk problem. We will begin in two space dimensions and find 
the Casimir energy of a Neumann line segment. For notational consistency with the 
three-dimensional case, we consider the x-z plane with the geodesic running along the 
^-axis and the segment running along the x-axis from x = —d to x = d. We define 
elliptical coordinates /x and 9 by 

x = d cosh ji cos 9 z = d sinh \i sin 9 (4) 

and then the boundary condition is applied at the points where the radial coordinate /i 
is zero. In elliptical coordinates, the problem is still separable, but there is no analog of 
angular momentum conservation. We still have separate angular and radial functions, 
but we now have the additional dimensionful parameter d. Thus the angular functions 
are now not simply functions of 9; they can now depend also on kd, where k is the wave 
number. Similarly the radial functions can depend on kd and r/d separately, rather 
than depending simply on the product kr. 

As described above, we have normal mode solutions that are even and odd under 
reflection of the z-axis. Since the line segment has Neumann boundary conditions, only 
the odd modes need to be modified from the free case. These become 

*I>oM, Q) = \ [e 2l5 Ho£(fi, q) + Ho%\fx, q)} se m (9, q) (5) 

where q = (dk/2) 2 . Here se m (9,q) is the odd angular Mathieu function and Ho$ and 
Ho$ are the corresponding odd radial functions of the third and fourth kinds, with 
Ho$ = Jo m + iYo m and Ho$ = Jo m — iYo m , where Jo m and Yo m are the radial 
functions of the first and second kind respectively. The phase shift 5(q) is defined by 
imposing the boundary condition, which yields 



HoS'(0,q) 



where the derivative is with respect to //. 

We have adopted the normalization conventions of [18], but have modified their 
notation to make the analogy to the circular case clearer. In these conventions the even 
and odd angular functions ce m (9,q) and se m (9,q) are normalized just like the cosm9 
and sinm^ solutions in the circular case, so that 

^ d9ce m (0, qf = 1^ d9se m (0, qf = vr . (7) 
Jo Jo 

As in the circular case, m = 0, 1, 2, 3 . . . for the even functions! and m = 1, 2, 3 . . . for 
the odd functions. The radial functions are then normalized so that they approach the 
corresponding Bessel functions at large radius. § Instead of singularities at zero radius 

J In the circular case, for to = the even solution is a constant, which is normalized to be -j= rather 

than cosO = 1, so that its normalization is consistent with the other modes. 

§ Note that the functions Je m and Jo m defined in [19] have an additional factor of 
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as in the circular case, the radial functions of the second kind have jump discontinuities 
— the singularity is now "spread" over the interfocal separation. 

From these wavefunctions, we form the normalized quantum field 



00 r°° k / - \ 

9) = £ / dk J- — (Je m (/i, q)ce m (0, q)bf + Jo m {n, q)se m (0, q)c^) e iu,t 

+ complex conjugate (8) 

where the odd term vanishes for m = 0. The stress-energy tensor for a minimally- 
coupled scalar field is 

T Xv = d A <R0 - l -r] Xv [d^d^] . (9) 

For a null vector, r)\ u V x V u = 0, so we have 

T Xv V x V v = {v a d a <t>f . (10) 

Taking our geodesic along the z-axis, V — (1, z), we have 

T Xv V x V v = 2 + {d^f . (11) 

We then substitute the result for the quantum field into this expression, subtract the free 
space result, and take the vacuum expectation value [20, 21, 22, 9]. For computational 
efficiency, we extend the k integration to the whole real axis and use contour integration 
to obtain an integral over the imaginary A;-axis, using k = in. We obtain [17] 

<0 2 > 4ff dK pf^ K 2 Ko m (^ V?se m (0, -V? • (12) 

^* m=l ^°m(0,y?) 

On the z-axis, terms with m even vanish, so we have 

1 00 r°° Tn' (0 in) 
<0 2 > = - E 7 *" J,in fW Konfr, y?) 2 S e m (vr/2, -v,) 2 (13) 

where the prime on the summation sign indicates that we sum over odd values of m. 
The sum over channels and integral over k are absolutely convergent, since we have 
subtracted the contribution of the free theory and we are away from the interactions, 
where all the remaining counterterms vanish. Similarly, on the z-axis we have 

7T 2 d 2 coshV i^i Jo K& n (Q,ipY 
In three dimensions, we employ oblate ellipsoidal coordinates, 



«^> 2 > - -^rai/r^izM^^^ 2 ^^ 2 '-^ 2 - (i4 » 



x = c^/(£ 2 + 1)(1 -^ 2 )cos0 y = d^ 2 + l)(l-r7 2 )sin0 z = (15) 

using the conventions of [23]. As in spherical coordinates, the solutions are indexed 
by n — 0, 1, 2, 3 . . . and m = — n . . . n. The angular solutions are spheroidal harmonics 
y™(ic; 77, 0), where c = fed, cos" 1 rj is the polar angle, and is the azimuthal angle. These 
are normalized analogously to ordinary spherical harmonics. The radial functions are 
R™^(ic; and R™( 2 \ic; where £ is the radial coordinate. These are normalized 
so that they approach the usual spherical Bessel functions at large radius. The factors 
of rki convert from prolate to oblate coordinates. 
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Now we consider Neumann boundary conditions on the disk £ = 0. If m + n is 
even, the functions obey the boundary conditions already. Otherwise we need to take a 
combination of iC (1) andiC (2) - With iC (3) = R™ {1) +iR" l{2) and iC (4) = R™ {1) -iR™ {4) 



we can write the desired radial function 



-<0 = ^ [e 2l5(lc) ^ (3) (<c; -«) + ^ (4) (<c; -ifl] (16) 

with the condition 

^ = iC (4)/ feQ) fl7) 

i^ (3) '(ic;0) 

where the derivative is with respect to the second argument. 

Passing to the imaginary A;-axis in the same way as in the two-dimensional case, 
we obtain [17] 

1 oo n f00 Dm(l)// n\ 

<* 2 > = -^£ E 7 d«« 3 % ; 7 '" fe (7;^^)l 2 iC (3) (7;-^) 2 (18) 

^ „=0m=-n J ° Rn '(75 0) 

where 7 = ic = i/cd = —red and the prime on the summation sign means that only odd 
values of m + n are included. On the axis, we have 



<0 2 > = -- E 7°° dK K * S| 7; °1 ^ ^)1 2 ^ (3) (7; -iO 2 (19) 

^ n =l J0 Rn ; '(7; 0) 

where we have specialized to m = because the contributions from nonzero m vanish 
on the axis, leaving only a sum over odd values of n. Similarly, on the axis we have 



<(<^) 2 > 



*4^l35(7;M)W(7;-ie)'(20) 

where the primes on the radial functions indicate derivatives with respect to the second 
argument. 

Once we have computed the result for a Neumann segment and disk, we must 
simply combine with the Dirichlet mirror result (see for example [24]), 

. in two dimensions, and 

<0 2 > + <(<9,0) 2 > = { 32 \ z3 . (21) 

in three dimensions 



1Qtt 2 z 4 



3. Results 



We carried out the two-dimensional calculations using the C++ Mathieu function 
package of Alhargan [25, 26], with minor enhancements to accommodate the extreme 
range of Mathieu functions needed to accurately compute this sum. We have also 
adapted the code to use our normalization conventions instead of those of [19]. The 
sums and integrals are then done by calling the C++ code from Mathematica routines. || 
In three dimensions, we used the Mathematica spheroidal packages of Falloon [27], with 

|| Mathematica does provide built-in support for both radial and angular Mathieu functions, but only 
for functions of the first kind (as of version 5.2). 
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Figure 4. Contributions to NEC in two dimensions (left) and three dimensions (right) 
for a Dirichlct plate with a hole of unit radius, as functions of distance along the axis 
passing through the center of the hole. Extrapolation is used for points at a distance 
less than 0.15 in the left panel and 0.25 in the right panel. The dotted lines show the 
perfect mirror result. 



minor enhancements to avoid memory leaks, allow compatibility with current versions 
of Mathematica, and improve efficiency for our application. Results in two and three 
dimensions are shown in Figure 4. We see that in both cases the hole has a dramatic 
effect, overwhelming the NEC- violating behavior away from the plate so that ANEC is 
obeyed. Although all field theory divergences are well under control, this calculation is 
still highly nontrivial for points near the hole, because both the perfect mirror and the 
disk have energies that diverge like l/z n+1 where n is the space dimension and z is the 
distance to the hole. The final result, however, is perfectly finite — the origin is just 
a point in empty space. As a result, we have to stop our calculation a fixed distance 
away from the origin, depending on our numerical precision (and patience). We then 
extrapolate the result to zero, and verify that it goes to a finite result with zero slope 
(without building this requirement into the extrapolation). Far away, the calculation 
approaches the perfect mirror result, which is also shown in Figure 4. 

Since the Dirichlet plate with a hole obeys ANEC, and the Neumann and Dirichlet 
cases typically contribute with opposite signs, we might then expect a Neumann plate 
to violate ANEC, with a negative contribution from the hole overwhelming the positive 
contribution from far away. It is straightforward to repeat this analysis for that case 
[17]. However, as shown in Figure 5, we see that ANEC is also obeyed for a Neumann 
plate with a hole. Furthermore, both the Dirichlet and Neumann results extend to the 
case of two plates in the limit of both large and small holes [17]. 

4. Conclusions 

Energy condition violation by Casimir systems initially looks dramatic, but in consistent 
field theory models it is modest or nonexistent. We do not have an example of 
ANEC violation in flat space — that is, due to quantum effects in a background of 
non-gravitational quantum fields. Violations are known to exist [28, 29, 30, 31] in 
curved space, when one considers the quantum effects of curvature caused by distant 
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Figure 5. Contributions to NEC in two dimensions (left) and three dimensions (right) 
for a Neumann plate with a hole of unit radius, as functions of distance along the axis 
passing through the center of the hole. Extrapolation is used for points at a distance 
less than 0.11 in the left panel and 0.25 in the right panel. The dotted lines show the 
perfect mirror result. 



masses. But the magnitude of this violation is typically much smaller than competing 
effects associated with the source of the curvature itself (unless we consider Planck- 
scale objects, where classical general relativity is unreliable), so it is not clear that 
this violation will persist in the full theory. If ANEC (possibly with appropriate 
modifications in curved space) is always obeyed by realistic quantum field theories with 
no uncontrolled external agents, it would prevent the appearance of exotic phenomena. 
Work is underway to investigate the possibility of finding analytic arguments to ensure 
that ANEC is always obeyed in flat space. 
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